CU-TP-806 
CAL-623 

IASSNS-AST-96/65 
astro-ph/9612209 

The Electron- Screening Correction for the Proton-Proton 

Reaction 

John N. Bahcall* 

School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 0854-0 

Xuelei Chenl^ and Marc Kamionkowski-'- 

Department of Physics, Columbia University, 538 West 120th Street, New York, NY 10027 

(December 1996) 

Abstract 
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solution agrees with the numerical solution to O(10~^). 

PACS number(s): 97.10.Cv,96.60.Kx,95.30.Cq 



Typeset using REVT^ 



*jnb@ias.edu 

^xuelei@phys. columbia.edu 
■i- kamion@phys.columbia.edu 



1 



I. INTRODUCTION 



Over the past three decades, much work has been devoted to refining the input data used 
in calculating solar-neutrino fiuxes. The comparison between the predicted and the observed 
fiuxes has important implications for particle physics and astrophysics. Most recently, a 
great deal of attention has been paid to a relatively minor effect, the electron screening of 
nuclear fusion reactions [|T]-|§] , and it has even been argued |0 that the discrepancy between 
observations and theoretical predictions might be reduced significantly if a different screening 
correction is adopted. 

We test in this paper the robustness of the standard WKB analytic treatment due to 
Salpeter []T| by solving numerically the relevant Schrodinger equation, including a Debye- 
Huckel screening potential, for the fundamental proton-proton (pp) reaction. The unscreened 
rate of this reaction can be calculated precisely using standard weak-interaction theory, 
accurate laboratory data for the two-proton system, and different refined deuteron wave 
functions in agreement with a variety of nuclear-physics measurements. Radiative correc- 
tions are also included in the standard calculation 

In the next Section, we re-derive Salpeter's analytic result for the weak-screening limit 
using a kinetic-theory approach (rather than Salpeter's thermodynamic arguments). In Sec- 
tion III, we calculate the proton-proton wave function in both a screened and unscreened 
Coulomb potential by numerical solution of the Schrodinger equation. We then use these 
results to evaluate numerically the electron-screening correction to the proton-proton reac- 
tion, thereby testing the validity of the standard WKB calculation. In an Appendix, we 
calculate a correction to Salpeter's result and find it to be negligibly small. 



II. REVIEW OF SCREENING CORRECTION 

To begin, we re-derive Salpeter's screening correction. To do so, we use kinetic theory 
to calculate reaction rates in both a screened and unscreened plasma. Although Salpeter's 
derivation was based on a thermodynamic argument, this alternative approach recovers the 
same results, and it will be useful for understanding the numerical work in the following 
Section. With our analytic approach, a very small correction to Salpeter's results is obtained 
and presented in the Appendix. 

The nuclear fusion rates in the solar interior are controlled primarily by Coulomb barriers. 
Therefore, the energy dependence of the fusion cross section is usually written as 

a{E) ^ ^exp(-2vrr/), (1) 
where S{E) is a function that varies smoothly in the absence of resonances, and 

Here, Zie and are the charges of the two colliding nuclei, and v is their relative velocity. 

The controlling factor, exp(— 27r?7), in Eq. ([1|) takes into account the probability for the 
nuclei to tunnel through the Coulomb barrier. It is obtained from the Coulomb potential 
V{r) = ZiZ2e^/r through the WKB approximation. 
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T{E) = exp £ mVcouiir) - E)Y'^dr 



exp(-2Z,Z.^^|/;^i-l du\ (3) 

= e-''^\ (4) 

where Tc is the classical turning point, defined by Vcoui(^c) = E, E is the kinetic energy, and 
fi is the reduced mass, 

17111712 ,_s 

fi = ■ . (5) 

mi + m2 

Here, mi and 1712 are the masses of the reacting nuclei. 

If the energies of the reacting nuclei have a Maxwell-Boltzmann distribution at a tem- 
perature T, the thermally-averaged cross section times relative velocity is 0] 



However, in the stellar interior each nucleus, even though completely ionized, attracts neigh- 
boring electrons and repels neighboring nuclei; thus, the potential between two colliding 
nuclei is no longer a pure Coulomb potential, but a screened potential Vsc{r). In the weak- 
screening case, the Coulomb interaction energy between a nucleus and its nearest few elec- 
trons and nuclei of the gas is small compared with the thermal energy ksT. In this case, 
the surrounding electrons and ions are only slightly displaced, and we obtain a screened 
potential of the form [|I| 

Vsc{r) = ^^^^^ exp(-r/rD). (7) 



Here, 

rn = C'^'[^-^j a, (8) 
is the Debye radius for the cloud; A = AiA2/{Ai + A2) is the reduced mass in atomic mass 



units; ( = \jY,i{XiZf / Ai + XiZi/Ai); Xj, Zi, and Ai are the mass fraction, charge, and mass 
number, respectively, of nucleus i; 

p-i/3(0.51 X lO^^cm) (9) 



is a measure of interparticle distance; p is the density in units of g cm~^; and Nq is Avogadro's 
constant. 

For the screened potential, the penetration factor is then given by 

/_2r I /•! rl 1^/^ \ 

r{E) = exp —^J2fiE / - exp (a;(l - u)) - 1 du] , (10) 
\ ri ^ Jo lu J / 
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where x = x{E) = rdro- Here, Tc is the classical turning-point radius defined by Vsc{T(^ = 
E. However, if x is small, then for the screened potential is roughly that for the unscreened 
potential: Tc — ZiZ2e^ / E. By expanding the exponential in the small-x limit (to be justified 
below), we obtain 

T{E) = exp[-2nr]{l - x/2)] = e'^'^^e^"^. (11) 

Although XTCT] depends on the energy, the effect of the correction on the thermally-averaged 
cross section can be approximated by evaluating xttt] at the most probable energy of inter- 
action, 

Eo = [{T^aZ^Z^kBTfiniAc^ /2)YI^ (12) 
= l.22QA{ZlZlATlYl^ keV (13) 

where m is the atomic mass unit, and Tg is the temperature in units of 10^ K. Then, 

{^v) J— ^^-^/o dES{E) exp(-27rry - E/ksT), (14) 

where the Salpeter factor /o is given by 

/o = e-o-^ = exp(0.188ZiZ2Cp'/'T6-'/'), (15) 

where p is the density in units of g cm~^, and 

xo = x{Eo) = 0.0133 {ZiZ2y/''A-'^''p'^^T^''^\. (16) 

For the pp reaction, xq ~ 0.01, which justifies the small-x approximation used above. Equa- 
tions (|6|)- (p!6D provide an alternative derivation of the Salpeter fl]] weak- screening formula. 



III. NUMERICAL RESULTS 

We now calculate numerically the cross section for the pp reaction for a Coulomb potential 
and a screened Coulomb potential to compare with the WKB calculation of the screening 



correction. To do so, we note that the reaction rate is proportional to A [|T0|, where A is 



the overlap integral of the proton-proton wave function and the deuteron wave function. 
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A = \ ^ / Ud{r)upp{r)dr, (17) 



where Op is pp scattering length, 7 = ^j2pE'i is the deuteron binding wave number, and E^ 
is the deuteron binding energy. The function Ud{r) is the radial part of the S'-state deuteron 
wave function. Our calculation in this Section follows the approach and notation of Ref . . 

For the purposes of this exercise, we use the McGee wave function |]lT| for the deuteron. 
If another wave function (which fits the deuteron data) is used, the overlap integral changes 
only slightly. Since we are here only investigating the effect of the screening correction to 
the reaction rate, our specific choice of the deuteron wave function is unimportant. 
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The radial wave function Upp{r) satisfies the radial Schrodinger equation, 

d^u r 1 



where 



u = -k^u, (18) 



R = ^ = 28.8198 fm, (19) 

k = nv/h is the center-of-mass momentum, and Vnuc('") is the short-range nuclear potential. 
For Vnuc{f) we use an exponential potential which yields the observed value for the scattering 
length and effective range [§]. Again, the overlap integral turns out to be practically inde- 
pendent of the detailed shape of the nuclear potential (as long as it matches the measured 
scattering length and effective range), so the choice of nuclear potential is unimportant for 
determining the screening correction. 

In the weak-screening case, Eq. (0) is replaced by 

' ' u = -k^u. (20) 



dr"^ 



Rr 



The solution to the Schrodinger equation is unique once the two boundary conditions are 
given. The first condition is u{0) = 0. The other boundary condition is obtained by noting 



that the asymptotic behavior of the wave function for r ^ rr, must be [12 



^p^V) ~ N^oni Sin (kr - ^ log(2A;r) + 5r') , (21) 
for the Coulomb potential, and 

u;i{r)^N,,sm{kr-Sl'), (22) 

for the screened potential, where 6'^°'^^ and 6^^ are phase shifts. Fixing the incident flux of 
protons for the Coulomb and the screened-Coulomb interactions requires A^^coui = A^sc- 

To solve this boundary-value problem, we integrate Eqs. (|18D and (^) from r = with 
the condition u{0) = and u'{0) = 1 to a large distance (about 10 r^), and then test that 
the solutions converge to the form of Eqs. (^) and (|2^) , respectively. From these numerical 
solutions, we obtain the normalizations N^^ and Ncoui- We then use the calculated wave 
functions to evaluate the overlap integral in Eq. ([T7|) both with and without screening. 
By squaring the ratio of the two overlap integrals, we determine numerically the screening 
correction to the cross section for the pp reaction. 

In Fig. we plot the WKB (dashed curve) and numerical (solid curve) results for the 
screening correction for the pp reaction as a function of relative momentum k. Our re- 
sults show that the discrepancy is small. For values of k at which the pp reaction occurs 
(k ^ 0.016 fm~^), the fractional difference is O(10~^). In fact, at smaller k, the fractional 
difference is expected to be even smaller, as argued below. The increased discrepancy at 
smaller k shown in Fig. |1| is due to numerical error in our calculation: accurate integration 
of the Schrodinger equation becomes increasingly difficult at smaller k since the asymptotic 
forms in Eqs. (PT|) and (521) are reached at progressively larger r. 
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FIG. 1. The screening enhancement factor as a function of k. The dashed curve is that for the 
WKB result, and the solid curve is that for the numerical result. 



Our main result is that a numerical solution of the screened Schrodinger equation for the 
proton-proton reaction reaction gives results in excellent agreement [to O(10~^)] with the 
rate calculated analytically using the usual WKB approximation, as originally formulated 
by Salpeter. 
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In this Appendix, we calculate a correction to Salpeter's screening formula and find it to 
be negligibly small. The integral in Eq. ( p!4D is usually evaluated by expanding in a power 
series of the inverse of a large quantity r, 



IV. CONCLUSION 
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APPENDIX: A SMALL CORRECTION 
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r = 3Eo/kBT = A2A87{ZfZ^AT^^y/\ 
The average product can then be written in a compact form P] 



{(Tv) = 1.3005 X 10" 



-15 



Z1Z2 

'at! 



1/3 



/oS'efrexp(-r) cm s 



3„-l 



where S^s is expressed in keV-barns. To first order in r ^ 

5 hS'EQ S"El 



5eff = S{Eo) 1 + 



12 



2S 



S 



E=Eo/ 



Expressing the various quantities in terms of their values at = 0, we find |T3 
S,^{E,) ~ S(0) 



1 + 



12r 



+ 



S'{Eo + ffcsT) , S"E^ fEo , 89 



S 



+ 



S 



2 72 ^ 



(23) 



(24) 



(25) 



(26) 



J E=0 



More accurately, however, we should include the factor e^'^^ in the thermal-average inte- 
gral, Eq. (|ll]). To do so, we rewrite the integral as 



(av) 



POO 

/ dES{E)exp{-2Trr] - E/keT + XTTT]). 
Jo 



I TXIl{kBTY JO 

Introducing the dimensionless quantity z = E/Eq, the exponential can be written as 



2r 

2717] — E/ksT + xnr] = z' 

3 



-1/2 + 

3 3 



(27) 



(28) 



To first order in Xq, the minimum point of the exponent is thus shifted to z = 1 — Xq, or 
E = Eq{1 — Xq). Using Laplace's method for asymptotic expansion of integrals, we find 
that the only 0(xo) correction to Eq. (|2^ ) is in the expression for S^s- To order 0{t^^,Xo), 
it is 



^eff = SiE) 1 + r 



-1 



^ 5S'{E)Eo , S"{E)Ei 
12 



2S{E) 



SiE) 



(29) 



^E=Eoil-x)^ 



In other words, to first order in xq, Eq in Eq. ( pS] ) should be replaced by £"0(1 — Xq). 
Expressed as S{0), we have 

^ ^ ^ _^ S'mEoil-x) + §kBT) _^ S"{0)Eo fE^^^ 



eff 



5(0) 



12r 



S{0) 



5(0) 



89 

2x) + ^ksT 



(30) 



where we have neglected terms of order 0{xq/t). We see that there is an 0{x) correction 
to the S' and S" terms. Since x is small (~ 10"^ for the pp reaction at the core of the 
Sun), and the S' and S" terms are generally small compared with the lowest-order term, 
these corrections are very small, typically ~ 0.1%. Therefore, the standard multiplicative 
correction factor (/o) should give a screened interaction rate which is accurate to 0(1%) in 
the weak-screening regime. Furthermore, since Xq increases only very slowly with increasing 
mass number, the standard correction should also be accurate for other fusion reactions 
which are in the weak- screening regime. 
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